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Abstract 

The charged hadron multipUcity fluctuations are considered in the canonical ensemble. The microscopic 
correlator method is extended to include three conserved charges: baryon number, electric charge and 
\ strangeness. The analytical formulae are presented that allow to include resonance decay contributions to 
■ correlations and fluctuations. We make the predictions for the scaled variances of negative, positive and all 
. charged hadrons in the most central Pb+Pb (Au+Au) collisions for different collision energies from SIS and 
^ AGS to SPS and RHIC. 
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I. INTRODUCTION 



The statistical models have been successfully used to describe the data on hadron multiplicities 



in relativistic nucleus-nucleus (A+A) colhsions (see, e.g., Ref. llEjS] and recent review |^). The 
applications of the statistical model to elementary reactions and/or to rare particles production 
have stimulated an investigation of the relations between different statistical ensembles. In A+A 
collisions one prefers to use the grand canonical ensemble (GCE) because it is the most convenient 
one from the technical point of view. The canonical ensemble (CE) P, ^, 0, 3 or even 

the microcanonical ensemble (MCE) ll| have been used in order to describe the pp, pp and e"^e~ 
collisions when a small number of secondary particles is produced. At these conditions the statistical 
systems are far away from the thermodynamic limit, the statistical ensembles are not equivalent, and 
the exact charge, or both energy and charge conservation, laws have to be taken into account. The 
CE suppression effects for particle multiplicities are well known in the statistical approach to hadron 



production, for example, the suppression in a production of strange hadrons j8|, antibaryons j^, and 
charmed hadrons 10] when the total numbers of these particles are small (smaller than or equal to 
1). The different statistical ensembles are not equivalent for small systems. When the system volume 
increases, V — > oo, the average quantities in the GCE, CE and MCE become equal to each other, 
i.e., all statistical ensembles are thermodynamically equivalent. 
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The fluctuations in high energy nuclear collisions (see, e.g., Refs. 
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25| and references therein) reveal new physical information and can be closely related 
to the phase transitions in the QCD matter. The particle number fluctuations for relativistic systems 
in the CE were calculated for the flrst time in Ref. |2Q] for the Boltzmann ideal gas with net charge 
equal to zero. These results were then extended to quantum statistics and non-zero net charge in the 



fsee also Refs. 
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CE |22l, |2a, 1221, I and to the MCE 



n 



and compared to the corresponding results in the GCE 



. l34|). Expressed in terms of the scaled variances, the particle number fluctuations 
have been found to be suppressed in the CE and MCE comparing to the GCE. This suppression 
survives in the limit V — oo, so the thermodynamical equivalence of all statistical ensembles refers 
to the average quantities, but is not applied to the scaled variances of particle number fluctuations. 

The aim of the present paper is to extend a microscopic correlator method to treat the hadron- 
resonance gas within the CE formulation. In Section |n] we calculate the microscopic correlators 
in relativistic quantum gas. This allows one to take into account Bose and Fermi effects, as well 
as an arbitrary number of the conserved charges in the CE. We also argue that the microscopic 



in the 



correlator approach gives the same results as the explicit saddle point CE calculations 
large volume limit V oo. In Section ITTTl we deflne the generating function to include the effects 
of resonance decays. This gives the analytical expressions for resonance decay contributions to the 
particle correlations and fluctuations within the CE and MCE. In Section HVl we calculate and make 



the predictions for the scaled variances of negatively, positively and all charged particles in central 
Pb+Pb (Au+Au) collisions along the chemical freeze-out line at different collision energies. Section 
El presents our summary and conclusions. Some details of the calculations are given in Appendix. 



II. EXACT CHARGE CONSERVATIONS IN STATISTICAL SYSTEMS 
A. CE Microscopic Correlator 

Let us consider the fluctuations in the ideal relativistic gas within the CE. Our primary interest 
is to include different types of hadrons, while keeping exactly fixed the global electric (Q), baryon 
(B), and strange (S) charges of the statistical system. The system of non-interacting Bose or Fermi 
particles of species i can be characterized by the occupation numbers j of single quantum states 
labelled by momenta p. The occupation numbers run over j = 0, 1 for fermions and n^^j = 
0, 1, 2, ... for bosons. The GCE average values and fluctuations of Up^i equal the following [s^: 

i^p^i) = r/ , ^ 1 ' (1) 



exp 



<i = = (K,.-K.))') = K,)(i + 7.K.)) • (2) 

In Eq. dU, T is the system temperature, mj is the mass of i-th particle species, 7^ corresponds to 
different statistics (+1 and —1 for Bose and Fermi, respectively, and 7^ = gives the Boltzmann 
approximation), and chemical potential /ij equals: 

= Qi f^Q + h fiB + Si fis , (3) 

where g^, bi, Si are the electric charge, baryon number and strangeness of particle of species i, 
respectively, while fig, fiB, fJ'S are the corresponding chemical potentials which regulate the average 
values of these global conserved charges in the GCE. 

The average number of particles of species i, the number of positive, negative, and all charged 
particles are equal: 

T A /"OO 

(iv.) ^ = / p'dp M , (4) 

i,qi>Q i,qi<0 i,<Ji^O 

where V is the system volume and gi is the degeneracy factor of particle of species i (a number of spin 
states). A sum of the momentum states is transformed into the momentum integral, which holds in 
the thermodynamic limit \^ — > 00. 



The microscopic correlator in the GCE reads: 



(6) 



where j is given by Eq. 0. This gives a possibihty to calculate the fluctuations of different 
observables in the GCE. Note that only particles of the same species, i = j, and on the same level, 
p = k, do correlate in the GCE. Thus, Eq. © is equivalent to Eq. only the Bose and Fermi 
effects for the fluctuations of identical particles on the same level are relevant in the GCE. 

In order to include the effect of exact conservation laws, we introduce the equilibrium probability 
distribution W{Anp^i) of the deviations of different sets {^p j} of the occupation numbers from their 
average value. In the GCE each Arip j fluctuates independently according approximately to the Gauss 
distribution law for Aup^i with mean square deviation (An^ J 



p,i 



Wg,c.e.{Anp^i) oc JJexp 

p,i 



(Anp^iY 

p,t 



(7) 



To justify Eq. (0) one can consider (see Ref. [l2|) the sum of rip^i in small momentum volume {Ap)^ 
with the center at p. At flxed (Ap)^ and V ^ oo the average number of particles inside (Ap)"^ 
becomes large. Each particle conflguration inside (Ap)^ consists of (Ap)^ ■ gV / {27i)^ ^ 1 statistically 
independent terms, each with average value (up^i) © and variance fpj (0). From the central limit 
theorem it follows then that the probability distribution for the fluctuations inside [Ap)^ should be 
Gaussian. In fact, we always convolve rip^i with some smooth function of p, so instead of writing the 
Gaussian distribution for the sum of Up^i in {ApY we can use it directly for Up^i. The next step is to 
impose exact conservation laws. The problem is to calculate the microscopic correlator with three 
conserved charges, Q, B, S, in the CE, i.e. when global charge conservation laws are imposed on each 
microscopic state of the system. The conserved charge, e.g., the electric charge Q, can be written in 
the form Q = YlpiQinp^i. An exact conservation law is introduced as the restriction on the sets of 
the occupation numbers {up^i}: only those sets which satisfy the condition AQ = J2pi(li'^^p,i = 
can be realized. Then the distribution should be modifled. This has been considered before for 
one conserved charge in the CE I27I and MCE . Now three charge conservation laws are imposed: 
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WceXAup^i) (X ]^exp 

dXqdXfydXs JJ^exp 
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It is convenient to generalize distribution 
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Pi* / 



2v^- 

p,i 



+ iXq QiArip^i + iXb hAup^i + iXs SiAn 



p,i 



using further the integration along imaginary axis in 



A-space. After completing squares one finds: 



WceX^rip/, Xg,Xb,\s) ocJJexp[ 



Pi* 



{Aup^i - Xgvliqj - Xbvlibj - KvliSj)' 



A A^ A^ 

+ y + y + Y^Ji^i - K^svliQi Si - XgXbvl^qi hi - XbXsvl J), Si 

The CE averaging takes the following form: 

^ ^ ^ r^ood\dXbdXs jT^Ylpidripi ... Wcej^^pX. Ag, Afe, A,) 
dXgdXbdXs Hp,* drip i W^.e. {AupX, Xg, Xb, A,) 

The CE microscopic correlator is as follows (see also Appendix): 



(9) 



(10) 



(Aup^iAnkj) 



c.e. ^P,i 



[qiQjMqg + bibjMbb + SiSjMss + {qiSj + qjSi) Mg^ - {qibj + qjbi) Mgb - {biSj + bjSi) Mi 
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where |A| is the determinant of the matrix, 

^A(g2) A{bq) A(sg)^ 
A= A{qb) A{b^) A{sb) , (12) 
^A(gs) A{bs) A(s2)^ 

with the following elements, A(g^) = J2pi1i'^p,i ? A{qb) = J2pi^i^i'^p,i ? ^^c. Mij are the corre- 
sponding minors of the matrix A, e.g.. 



det 



A{qb) A(62 
A{qs) A{bs) 



(13) 



In the case of conservation of only one (electric) charge, this reduces to |y4| = A(g^), Mgg = 1. To 
make these formulae more transparent we write one of the minors explicitely, 

M.. = A(g^) . A{b') - [A{qb)r = (j2q^ v^] ■ 6? vl) - (j^^A vl) , (14) 

\ P,i / \ k,j / \ p,i 



The sum, Y2pi 5 means integration over momentum p, and summation over hadron-resonance 
species i. The microscopic correlator can be also used in the MCE. The exact energy conserva- 
tion is imposed with 5 (AE) = 5 (^^p j a/ m? + Arip j j . This would lead to additional terms in the 
r.h.s. of Eq. dTT)) proportional to Xlp.il^^i "^p.i 1 Ylip,i a/ mf + qi Vp^i , etc. 

The microscopic correlator (fTT|) can be used to calculate correlations and fluctuations of different 
physical quantities in the CE. The first term in the r.h.s. of Eq. (jllj) corresponds to the microscopic 
correlator ^ in the GCE. The additional terms reflect the (anti) correlations among different parti- 
cles, i ^ j, and different levels, p ^ k, that appeared due to the global CE charge conservations. Let 



us concentrate on the particle number fluctuations. One can calculate the correlations in the GCE 
and CE, respectively, 

{ANi AN, ) = (Arip,, Aukj) = ^ v^^ , (AA^^ AA^, = ^ {Arip^i Ankj)c.e. • (15) 

p,k p p,k 

The CE scaled variance reads: 



<c.e. [ 1 



(iV.)e.e. 

- ^^JJ^ il'^Mgg + b'^Mbb + s'^.Mss + 2qiSiMgs - 2qibiMgb - 2hsiMbs) ] ■ (16) 

In Eq. flT^ we used the fact that {Ni)c.e. is equal to {Ni) ^ in the GCE aiV ^ oo, and introduced 
the scaled variance in the GCE, 

Note that the CE result ()16|) is obtained in the thermodynamic limit, and it does not include a 
dependence on V . Thus the method can not be used to obtain the finite volume corrections. A 
nice feature of the microscopic correlator method is the fact that particle number fluctuations and 
correlations in the CE, being different from those in the GCE, are presented in terms of quantities 
calculated within the GCE. 



B. Saddle Point Expansion Technique 



In this subsection we discuss the method of treating the CE at finite volume V . Let us for 
simplicity consider the CE with only one conserved charge, Q, and only one sort of particles with 
charges +1 and —1. The microscopic correlator method (jl6ll7|) then gives: 
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1 - 



±2 



and for zero value of the total net charge, Q = 0, this reduces to 



± 



L^p p 



±2 



2 ^g.ce. 



(19) 



Let us start with an example of Boltzmann approximation, 7 — 0. For neutral system, Q = 
one finds in the GCE: 



Zg.c.e. = exp (2z) 



{N±) 



{Nl) 



z + 



ijj 



g.c.e. 
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and in the CE aX V ^ 00 



3: 



Io{2z) , (Ar±), 



'-Tz 



(20) 
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where z = gV(27i'^) ^Tm?K2{m/T) is one particle partition function. It then follows: 



2 

C.e. r-u 



\ [1 + 0{llV)] - \ 



[22] 



also 



which coincides with Eq. fll9|) . This result has been obtained for the Boltzmann gas. To justify it 
for the Bose and Fermi gases we consider now a systematic sad dle-poin t expansion 0, (see 
3^, Is^). The CE partition function is defined as follows P 



ZceiQ) = E ) - 



N_) 



{n+,rap } 



2ti 



exp{-iQ(f)) Zg,c.e.{<P) , (23) 



with 



^g.c.e 



(0) = exp 



(- S^r^^ " 7A.e-^^/--^) + ln(l - ^A.e-^^/-'^)]) 



(24) 



/o 7 

where is the degeneracy factor, Sp = {p^ + m^)^/^, and 7 = +1 and —1 correspond to Bose and 
Fermi statistics, respectively, while the limit 7 — gives the Boltzmann approximation. The A+ and 
A_ in Eq. are auxiliary parameters that are set to one in the final formulae. A substitution of 
exp(±2(/)) in Eq. (j24|) by exp (±/i/T) leads to well known expression of the GCE partition function, 
Zg.c.e.-i with chemical potential j35|. In Eq. E = J2p^pi''^p + '^p) N± = X^p'^p- 
expands the logarithm in Eq. in the Taylor series, — 7"^ ln(l — 7 x) = Yl'^=il"'~^ x"'/n. This 
leads to: 



ZcAQ) 



2tt 



27r 



exp -tQ(p + 



9V_ 
27r2 
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2tt 



2tt 
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gv_ 

27r2 



1 °° i-i 



-Isp/T 



[(A+ e*<^)' + (A_ e-'I^Y] (25) 



ri=0 n=0 

The Boltzmann approximation, 7^0, corresponds to only one term, Z = 1, in the sum from Eq. ()25|). 
Using the following notations. 



-n = p'dpY. ^"-^"^ ^ 

-^0 Z=l «=1 



AL7 



27r2 / 



K2{lm/T) , 



(26) 
(27) 



where oc are the so called cumulants, one can easily get the following formula: 



ZceXQ) = exp(/t+ + Kq) / —exp -iQ(f) + -r ( « 



-(^0)" + «:;(-z0)") 



:281 



At A± = exp(±/x/T) the cumulants nf" give the GCE values: 



(29) 



The average values and fluctuations in the CE can be obtained as the following: 



A±=l 



A ^ 
dX± V ^ dX± 



A±=l 



To calculate (jHUj) one needs to estimate the following integrals, 



360^ 

At oo an integrand in (f3T|l has a strong maximum at = 0, which leads to the result: 

1 K4 



4^2 



The Eq. (001) then leads to 
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n=l 
00 

E 

n=l 



Zq 



i=l n=l 



At Q = 0, Eqs. dSSl) and (jSl read: 

(iV±)c.e. 



/€l 1 



/ 2 \ 

E-0-£r + f(v^-)) 

n=l ^ ^ 
ra=l \ 2 ^1=1 n=l ^ 



4Ki/t2 

00 00 / / r , \ 2 



(30) 
(31) 
(32) 

(33) 
(34) 

(35) 



The scaled variance equals: 



K2 
2k'i 



1 



(36) 



(37) 



which coincides with Eq. (fT^ in the large volume limit V ^ 00. 

One again observes that the global conservation laws lead to the correction to average particle 
numbers, (iV±)c.e. = {N±)[l — 0{1/V)]. It equals K3/(4kiK2) and leads to additional terms to 
(A^±)c.e. {N±)c.e. proportional to V. These terms, however, are cancelled out in the variances, 
and Eq. ^T^i obtained from the microscopic correlator remains valid. This gives justification of the 
microscopic correlator approach, which assumes the equality {np^i)c.e. = (np^i) in the thermodynamic 
limit — 00. 

III. EFFECT OF RESONANCE DECAYS 



A. Generating Function 



Resonance decay has a probabilistic character. This itself causes the particle number fluctuations 
in the final state. The average number of final particles from resonance decays, and all higher 



moments including particle correlations can be found from the following generating function: 

R \ r i J 

where 6^ is the branching ratio of the r-th branch, nf^ is the number of i-th particles produced 
in that decay mode, and r runs over all branches with the requirement Ylr^r^ ~ Note that 
different branches are defined in a way that final states with only stable (with respect to strong and 
electromagnetic decays) hadrons are counted. The Aj in Eq. (jHHj) are auxiliary parameters that are 
set to one in the final formulae. The averages from resonance decays can be found as the following: 

R ^ R r R 

iVWV^ . ^(N. iV,), ^ A. A (a, A g) 

= 5^ [ iV^ (iV^ - 1) {n,)R {n,)R + Njt {n, n,)n ] , (40) 

R 

where {rii nj)^ = J2r ^r'^i!r^f,r ■ The averaging, (■ ■ ■ )r, in Eq. (jH^ means the averaging over reso- 
nance decays. The formula (jHH|l originates from the fact that the normalized probability distribution, 
P{N^), for the decay of Nr resonances is the following: 

PiN'R) = Nr\ (E^^ - ' (41) 

where A'^Jj correspond to the numbers of -R-th resonances decaying via r-th branch. 
The scaled variance uj'^r due to decays of -R-th resonances reads: 



2 



To illustrate Eq. ()42|) some examples are appropriate. It follows from Eq. that cj}^ = if n^.^ 
were the same in all decay channels. The uj^^ also vanishes if there was only one decay channel, i.e. 
6f = 1. Let there be an arbitrary number of x-th type decay channels with nf^ = 1 and y-th type 
ones with n^y = 0. From Eq. (ji^ one finds tu)^ = 1 — 6^ > 0, where 6^ is the total probability of 
x-th type decay channels. If n^^ = 2 and n^y = 0, then c<j}^ = 2(1 — b^). In general, Eq. tells 
that resonance decays generate fluctuations of z-th hadron multiplicity if raf^ are different in different 
decay channels. If nf^ is larger than 1 in some of these channels, the fluctuations become stronger. 

The Eqs. (|39l4Up assume some fixed values of A^^^;. In a real situation, Nr fluctuate, and this is 
an additional source of the particle number fluctuations. One finds: 

^R = //.r\ \ = ^R + {ni)RUjR, (43) 



where resonances act as sources of particles, similar to the so called independent source model 
and the scaled variance, 



= ^ ""'^ ,^ , (44) 

corresponds to the thermal (GCE or CE) fluctuation of the number of resonances. 



B. Grand Canonical Ensemble 

The average number of i-particles in the presence of primary particles A^* and different resonance 
types R is the following: 

im = (iV;) + E(^«)E^X. = {Nn + J2^Nn){n,)n (45) 

R r R 

The summation YIr runs over all types of resonances. The (. . .) and (. . correspond to the GCE 
averaging, and that over resonance decay channels. 

From Eqs. (jHHl IH^ one finds the GCE correlators : 

{AN, AN,) = {AN* AN*) + [(AiV^^ {n,)R{nj)n + {Nr) {An,An,)n\ • (46) 

R 

The terms proportional to {NR){Njii) and {N*){Nr) cancel each other in the GCE calculations of 
the correlator PU)) . 



C. Canonical Ensemble 

All primary particles and resonances become to correlate in the presence of exact charge conser- 
vation laws. Thus for the CE correlators we obtain a new result: 

{AN, AN,),,,, = {AN* AN*),,,, + ^^{Nr) {An.An,)^ + ^(AiV* AiV^),.,. 

R R 

+ J2{AN* ANr),,,, {n,)R + J2{ANr ANr,),,,, {n,)^ {n,)^, . (47) 

R R,R' 

Additional terms in Eq. (jT7|) compared to Eq. (jiUI) are due to the correlations induced by exact 
charge conservations in the CE. The Eq. (j47p remains valid in the MCE too with (. . .),,,, replaced 

by (. . .)m.c.e.- 



IV. SCALED VARIANCES ALONG THE CHEMICAL FREEZE-OUT LINE 



In this section we present calculations of the CE and GCE fluctuations along the chemical freeze- 
out line in central Pb+Pb (Au+Au) for both primordial and final state distributions. 



At chemical freeze-out the hadronic gas is usually described by the following parameters: tem- 
perature T, chemical potentials {fJ^B, f^s^ fJ'Q)^ and strangeness suppression factor 75 to account for a 
undersaturation of the strange sector. The GCE has proven to be sufficient for thermal model anal- 
ysis of mean multiplicities in central Pb-Pb and Au-Au collisions at most colliding energies. Only at 
lower energies, where only a few strange particles are produced, the CE effects of an exact strangeness 
conservation become visible. This leads to the CE suppression of yields of strange particles when 
compared to the GCE. In the energy range discussed below this is only the case for the SIS data 
point. On the other hand, for multiplicity fluctuations the exact conservation laws are important for 
all colliding energies. 

Thermal model analysis has provided a systematic evolution of the parameter set with beam 
energy and size of colliding system and allows for phenomenological parametrization, giving the 
thermal model almost predictive qualities. A recent discussion of system size and energy dependence 



of freeze-out parameters and comparison of freeze-out criteria can be found in Refs. 



There are several programs designed for the statistical analysis of particle production in relativis- 
tic heavy- ion collisions, see e.g., SHARE j2| and THERMUS 37 1. In this paper an extended version 
of the THERMUS thermal model framework is used. With increasing colliding energy, the tem- 
perature increases and more energy for particle production becomes available. This is accompanied 
by a drop in fj,B, which can be parameterized by the following function 

, , , 1.308 GeV , 

= 1 + 0.273 ^ ' (''^ 

where the cm. nucleon-nucleon collision energy, ^/SNN, is taken in GeV units in Eq. PH|) . 

The electrical chemical potential /ig can be further adjusted to give the charge to baryon ratio of 
heavy nuclei, Q/B ^ 0.4. Strange chemical potential fis is constrained by requiring the system to 
be net strangeness free, S = 0. Finally the temperature is chosen to match a condition, (E)/ (N) ^ 
1 GeV 0|, for energy per hadron. In order to remove the remaining free parameter, 7^, we use the 
following parametrization |^: 

/ X 23 T\ 

75 = 1 - 0.396 exp '■ . (49) 



Numerical fitting functions allow to meet all the above criteria simultaneously and thus to choose a 
parameter set, (T, /i^), for each given collision energy. The corresponding chemical freeze-out line 
in the T — fiB plane is shown in Fig. 1. There is obviously some degree of freedom as to choose 
a particular parametrization for some parameter or value of the average energy per particle. This 
particular choice is in good agreement with thermal model fits done in Ref. |3|. The center of mass 
nucleon-nucleon energies, ^/SNN, quoted in Table I correspond to beam energies at SIS (2 AGeV), 
AGS (11.6 AGeV), SPS (20, 30, 40, 80, and 158 AGeV), and two top colliding energies at RHIC 
(^/^ = 130 GeV and 200 GeV). 



y/SNN [GeV] 


T[MeV] 


f^B[MeV] 


75 


PB [fm 3] 


2.32 


64.9 


800.8 


0.642 


0.061 


4.86 


118.5 


562.2 


0.694 


0.111 


6.27 


130.7 


482.4 


0.716 


0.117 


7.62 


138.3 


424.6 


0.735 


0.117 


8.77 


142.9 


385.4 


0.749 


0.115 


12.3 


151.5 


300.1 


0.787 


0.104 


17.3 


157.0 


228.6 


0.830 


0.088 


130 


163.6 


35.8 


0.999 


0.016 


200 


163.7 


23.5 


1 


0.010 



TABLE I: Chemical freeze-out parameters for central Pb+Pb (Au+Au) collisions. 




FIG. 1: Chemical freeze-out curve for central Pb+Pb (Au+Au) collisions. 



Once a suitable parameter set is determined, mean occupation numbers and fluctuations can be 
calculated using Eqs. dH) and The scaled variances of negative, positive, and all charged particles 
read: 



where 

HAN.f) = Yl i^NAN,) , ((AiV+)') = Yl i^NAN,) , (51) 

qi<0,qj<0 qi>0,qj>0 

{(AN^hf) = Yl i^NANj) . (52) 

For the primordial hadrons, mean multiphcities, (iVj), in Eq. ()50|) are given by Eq.(@]), and correlators, 
(ANiANj) , in Eqs. ()51I52|) are either the GCE or CE correlators from Eq. (fT3j) . For final state mean 
multiplicities, Eq. (j45|l is used, and correlators are calculated with Eq. (j46|l in the GCE or Eq. 
in the CE, respectively. For final hadrons the summation needs to be extended to all stable particles 
with corresponding charges and all unstable resonances which have these charged particles in their 
decay channels. Figures |21 El and HI show scaled variances for negatively charged particles, u~ , 
positively charged particles, u^, and all charged particles, u'^'^, respectively, as functions of y^^iViv- 
Four cases are considered, namely, primordial and final state particles in both the GCE and CE. The 
relevant primordial and final state values for various colliding energies are summarized in Tables ITTl 
and mil respectively. 
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L0 + 




^SNN [OeV] 


GCE 


CE 


GCE 


CE 


GCE 


CE 


2.32 


0.982 


0.373 


0.976 


0.115 


1.027 


0.775 


4.86 


1.027 


0.677 


1.013 


0.261 


1.059 


0.575 


6.27 


1.036 


0.737 


1.021 


0.296 


1.062 


0.564 


7.62 


1.041 


0.779 


1.027 


0.321 


1.065 


0.560 


8.77 


1.044 


0.808 


1.030 


0.339 


1.066 


0.558 


12.3 


1.049 


0.872 


1.037 


0.380 


1.066 


0.557 


17.3 


1.052 


0.929 


1.042 


0.418 


1.065 


0.559 


130 


1.054 


1.050 


1.053 


0.523 


1.056 


0.548 


200 


1.055 


1.053 


1.053 


0.529 


1.056 


0.545 



TABLE II: Primordial scaled variances in the GCE and CE for central Pb+Pb (Au+Au) collisions. 



The column pp in Table I allows for a comparison with previously reported values of primordial 



scaled variances 



29| (a good agreement is found). The standard THERMUS particle table includes 



all strange and light fiavored particles and resonances up to about 2.6 GeV. Only strong and elec- 
tromagnetic decays are considered, weakly decaying channels are omitted. It should be mentioned 
that, in particular, heavy resonances do not always have well established decay channels, thus there 
are always some ambiguities in the implementation of resonance decays in respective thermal model 
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^SNN [OeV] 
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GCE 


CE 


2.32 


1.048 


0.403 


1.020 


0.116 


1.025 


0.777 


4.86 


1.354 


0.848 


1.195 


0.327 


1.058 


0.621 


6.27 


1.421 


0.967 


1.201 


0.395 


1.068 


0.643 


7.62 


1.464 


1.059 


1.198 


0.449 


1.076 


0.668 


8.77 


1.491 


1.124 


1.194 


0.486 


1.082 


0.687 


12.3 


1.542 


1.268 


1.182 


0.571 


1.095 


0.732 


17.3 


1.576 


1.387 


1.171 


0.643 


1.105 


0.770 


130 


1.619 


1.613 


1.138 


0.802 


1.128 


0.824 


200 


1.620 


1.617 


1.136 


0.808 


1.130 


0.822 



TABLE III: Final state scaled variances in the GCE and CE for central Pb+Pb (Au+Au) collisions. 
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FIG. 2: The scaled variances for negatively charged particles, uj~ , along the chemical freeze-out line for 
central Pb+Pb (Au+Au) collisions (see Fig. Different lines present primordial and final GCE and CE 
results. 



codes. Details about the THERMUS decay convention can be found in Ref . [33| • The quoted values 
of the scaled variances are valid in the thermodynamic limit and assume that all charge carriers are 
detected. For high fiB (low collision energies) the multiplicity of positively charged particles, A^+, is 
enhanced in a comparison with A^„, while the fluctuations uj~^ are suppressed in a comparison with 
uj^. At vanishing net charge density (high collision energies), uj~^ and uj~ have the same asymptotic 
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FIG. 3: The scaled variances for positively charged particles, w^, along the chemical freeze-out line for 
central Pb+Pb (Au+Au) collisions (see Fig. Different lines present primordial and final GCE and CE 
results. 
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FIG. 4: The scaled variances for all charged particles, uj , along the chemical freeze-out line for central 
Pb-|-Pb (Au+Au) collisions (see Fig. Different lines present primordial and final GCE and CE results. 

values. The scaled variance for all charged particles u'^'^ has the same value in the GCE and CE for 
a neutral system, for both primordial and final state. The effect of resonance decays remains small 
at low collision energies (i.e. small temperatures), while becoming sizeable even at the lowest SPS 
energy. 



Some important qualitative effects are seen in Figs. El El and |31 The effect of Bose and Fermi 
statistics can be seen in primordial values in the GCE. At low temperatures Fermi statistics dominate, 
^t-c-e.^^t^c.e. < 1; while at high temperature (low hb) Bose statistics dominate, i^^c.e.; '^gx.e. > 1- 
the chemical freeze-out line, uj~^^ is always slightly larger than 1, as vr" is the dominant negative 
particle at low temperature too. A bump at small collision energies in ^^ce final particles is 
due to the decay into 2 positively charged hadrons, p + 7r+. This single resonance contribution 
dominates at small collision energies (temperatures), but becomes relatively unimportant at high 
collision energies. A minimum in uj~^ for final particles is seen in Fig. |21 This happens as a result 
of the following effects. Since the number of negative particles is relatively small, (A^-^) <^ at 
low collision energies, the CE suppression effects are also small. Low collision energies correspond 
to small temperatures of the hadron-resonance system, and the resonance decay effects are small 
too. With increasing ^Jsnn, the CE effects increase and this makes uj~^ smaller, but resonance decay 
effects increase too and they work in an opposite direction making uj~^ larger. A combination of 
these two effects, CE suppression and resonance enhancement, leads to a minimum structure of uj~^ 
seen in Fig. El 

The results for scaled variances presented in Figs. 2-4 and Tables II, III correspond to an ideal 
situation when all final hadrons are accepted by the detector. To compare our calculations to 
experimentally obtained values of uj the acceptance and resolution need to be taken into account. 
Observing only a fraction q of final state particles dilutes the effect of global charge conservation. 
Even though the primordial particles at the chemical freeze-out line are only weakly correlated in the 
momentum space, this is no longer valid for final state particles as the decay products of resonances 
are not re-thermalized. Neglecting the momentum correlations due to resonance decays (this is 
approximately valid for uj^ and uj~ , and much worse for uj'^^) the following a ppr oximation for the 
scaled variances of experimentally accepted particles can be used (see e.g., 20,|2a]), 

^acc = 1 - g + g t^47r , (53) 

where uj^t^ refers to an ideal detector with full 47r-acceptance. In the limit of a very 'bad' (or 'small') 
detector, g — 0, all scaled variances approach linearly to 1, i.e., this would lead to the Piossonian 
distributions for detected particles. 



V. SUMMARY 



The multiplicity fluctuations of hadrons in relativistic nucleus-nucleus collisions have been consid- 
ered in the statistical model within the canonical ensemble formulation. The microscopic correlator 
method, previously used for one conserved charge, has been extended to include three conserved 
charges - baryon number, electric charge, and strangeness. The analytical formulae for the resonance 



decay contributions to the correlations and fluctuations have been found. Using the full hadron- 
resonance spectrum we have calculated the scaled variances of negative, positive and all charged 
particles for primordial and final hadrons at the chemical freeze-out in central Pb+Pb (Au+Au) 
collisions for different collision energies from SIS and AGS to SPS and RHIC. Both the CE and 
resonance decay effects for the multiplicity fluctuations have been discussed. A comparison with the 
NA49 data in Pb+Pb collisions at the SPS energies can be done for the sample of most central events 
with the number of projectile participants being close to its maximal value, ■ /V|r °"' ~ A, to avoid the 
fluctuations of the number of nucleon participants (see discussion in Ref. j2J,|25[). These NA49 data 
will be available soon 3^. The predictions of the statistical model within the CE formulations can 
be done for uj~ and uj~^. In this case the experimental acceptance can be approximately introduced 
by a simple procedure based on Eq. We find a qualitative difference between the CE results, 

< 1, and the GCE ones, cu^c.e. > 1, for the accepted particles. 
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APPENDIX A 



The n-dimensional Gauss integral equals the following 



/+00 r~ 
-oo J — ( 



exp 



i,k=l 



dxi . . . dXr, 



where 



\A\ = detA = ^(-l)^+M,,fcM,,fc 



TT 



n/2 



(Al) 



(A2) 



1=1 

and Mj jfc = (— l)*+'^9|y4|/9^i ^ is a complementary minor of the element fc. One also finds: 
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